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摘 要 : 基于 2000 一 2014 年 的 MODIS 数据 和 Landsat 数据 ,利用 温度 植被 干旱 指数 反映 地 表 湿 润 状 况 , 分 析 了 退 


耕 还 林 


- 程 以 来 黄土 高 原 地 表 湿 润 状况 与 植被 覆盖 的 时 空 变 化 特征 及 其 相互 关系 。GD 2001—2014 年 黄土 高 原 归 
一 化 植被 指数 增 速 为 6.3% + (10a) ~ ,温度 植被 干旱 指数 降 速 为 -4.5% - 


(10a) ,植被 覆盖 与 地 表 湿润 状况 均 


旦 增加 趋势 ,分 布 呈 “东南 高 .西北 低 ” 的 特征 ,大 体 上 可 以 看 出 两 者 具有 正 的 空间 相关 性 。@@ 从 TVDI Al NDVI 的 
变化 趋势 来 看 ,黄土 高 原 植被 覆盖 与 地 表 湿 润 状况 整体 上 呈现 出 增加 趋势 。 空 间 分 布 上 反映 出 ,植被 覆盖 变化 基 


本 上 与 地 表 湿润 状况 变化 分 布 一 致 ,但 是 仍然 存在 区 域 性 差异 。(8) 从 相关 系数 来 看 ,黄土 高 原 植被 覆盖 与 地 表 湿 
润 状况 呈现 出 正 相 关 。 从 变化 趋势 对 应 关系 来 看 ,14 a 间 黄 土 高 原 东北 一 西南 一 线 山 地 地 区 植被 覆盖 与 地 表 湿 润 


状况 变化 有 利于 生态 环境 改善 。 人 类 活动 强烈 的 平原 地 带 极 不 利于 生态 环境 改善 。 值 得 注意 的 是 , 退 才 
区 是 未 来 生态 环境 变化 监测 和 管理 的 重点 区 域 ,以 免 引起 


还 林 重 点 


次 退化 。 
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植被 覆盖 变化 是 区 域 生 态 系统 环境 变化 的 重要 
指示 ,是 衡量 地 表 植 被 状况 的 重要 指标 ,对 水 文 . 生 
AS .全 球 变化 等 都 具有 重要 意义 ;同时 ,植被 覆盖 又 
是 影响 土壤 侵蚀 与 水 土 流失 的 主要 因子 。 因 此 ,从 
保护 和 改善 西部 生态 脆弱 区 生态 环境 出 发 实施 退耕 
还 林 草 政策 ,将 水 土 流失 严重 的 坡 耕 地 和 土地 沙化 
严重 的 耕地 停止 耕种 ,造林 种 草 ,建设 林 草 植被 。 自 
1999 年 退耕 还 林 草 政策 实施 以 来 ,黄土 高 原 迅 速 扩 
大 人 工 林 建设 ，1999 一 2010 年 累计 造林 面积 
1.89 x107 hm?!) ,植被 覆盖 率 从 1999 年 的 31.6% 
增加 到 2013 年 的 59.6%" ,增幅 达 88% ,植被 覆盖 
SEA BI IAS? ;入 黄河 泥 沙 量 锐 减 ,2013 年 
已 达到 了 与 唐 代 相 当 的 水 平 (2 Gt . a” ) ,并 趋 于 平 
Raa) ,黄土 高 原水 土 流失 得 到 缓解 ,生态 环境 
不 断 改善 ,鉴于 此 ,国家 "十 三 五 规划 建议 扩大 退 
耕 还 林 草 政策 的 实施 。 但 是 , 相 比 于 天 然 植被 ,人 工 
植被 光合 作用 消耗 的 水 分 较 大 ,加 之 前 期 退耕 还 林 
草 植物 生长 所 需 水 分 的 增加 ,引起 干旱 一 半 干 旱地 
区 地 表 湿 润 状况 降低 ,土壤 容易 出 现 早 化 特征 ,增加 
了 植被 衰败 的 可 能 性 。 随 着 气温 的 不 断 升 高 ,地 
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表 蒸 散发 加 大 、 河 湖 萎缩 与 干 润 , 气 候 和 环境 的 干旱 
化 成 为 严峻 的 生存 环境 问题 ,进而 使 得 地 表 湿 润 状 
况 继续 减 小 ,土壤 时 化 程度 加 重 , 不 仅 反 馈 给 人 工 植 
被 发 生 水 分 胁迫 ,致使 人 工 植 被 正常 生理 代谢 失调 ， 
从 而 影响 地 表 植 被 覆盖 的 变化 ,而 且 当 土壤 形成 不 
可 道 转 的 干 层 时 ,又 会 影响 到 降水 或 地 表 水 入 渗 补 
给 地 下 水 ,以 及 陆地 水 循环 *”” 。 如 果 和 忽略 了 生态 
脆弱 区 本 底 自 然 环境 特征 及 土壤 水 分 时 空 分 布 特 
征 , 育 目 扩大 生态 恢复 和 植被 建设 面积 ,可 能 引起 人 
工 植被 退化 ,其 至 导致 群落 衰败 和 生态 系统 退化 , 造 
成 巨大 生态 经 济 损失 ,限制 了 生态 恢复 项 目的 可 持 
续 发 展 ,与 生态 工程 预期 成 效 产生 差距 。 因 此 ,监测 
退耕 还 林 工 程 以 来 地 表 湿 润 状 况 与 植被 覆盖 变化 特 
征 ,以 及 分 析 两 者 的 相互 关系 ,进一步 探究 生态 工程 
的 实施 是 否 利 于 生态 环境 改善 ,具有 重要 的 理论 和 
现实 意义 。 

地 表 湿 润 状况 的 传统 监测 方法 有 效 性 较 高 ,是 
地 表 湿 润 状况 监测 .研究 、 应 用 的 主要 手段 ,但 是 该 
方法 采样 速度 慢 .耗费 大 量 人 力 和 物力 、 使 用 范围 有 
限 ,加 之 测 站 分 布 不 均 ,难以 实现 大 范围 实时、 精确 
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的 监测 。 随 着 遥感 技术 的 迅速 发 展 , 遥感 数据 的 多 
时 相 、 多 光谱 多 分 辩 率 ,高 光谱 等 特性 ,使 得 地 表 湿 
润 状况 快速 .实时 OAS WE TM OW AT AB" 。 在 地 
表 湿 润 状况 遥感 反 演 中 ,通过 地 表 温 度 和 植被 覆盖 
特征 空间 反映 地 表 湿 润 状况 的 遥感 监测 得 到 了 广泛 
应 用 "~"。 由 于 温度 植被 干旱 指数 (temperature 
vegetation dryness index, TVDI) 计算 过 程 输入 数据 
> 简单 易 行 .生物 物理 学 含义 明确 等 特点 ,因此 本 
研究 利用 植被 覆盖 与 地 表 温 度数 据 构 建 了 TVDI, 用 
以 反映 退耕 还 林 以 来 黄土 高 原 地 表 湿 润 状况 ,以 期 
为 后 续 退 耕 还 林 工 程 的 合理 实施 提供 一 定 的 依据 。 


1 资料 与 方法 


1.1 资料 来 源 

本 文 所 用 数据 包括 黄土 高 原 2001—2014 年 
MODIS-NDVI( normalized difference vegetation index, 
NDVI) .MODIS-LST (land surface temperature, LST) 
LAR 2001 年 和 2014 年 Landsat 数据 。MODIS-NDVI 
来 自 LP DAAC 中 的 MODIS 产品 的 MOD13A2 数据 
集 ,空间 分 辨 率 为 1000 m, 时 间 分 辨 率 为 16 d ,利用 
单位 转换 系数 0.000 1 将 其 转 为 0 ~1, 采 用 最 大 值 
合成 法 将 NDVI 数据 转换 为 年 际 数据 ,以 反映 逐年 
植被 覆盖 状况 ;MODIS-LST 来 源 于 MOD11A2 ,空间 
分 辨 率 为 1 000 m, 时 间 分 辩 率 为 8 d, 利 用 单位 转换 
系数 0. 02 将 其 转 为 开尔文 单位 (Kelvin) ,由 于 LST 
存在 大 面积 的 无 效 数 据 , 所 以 采用 Cressman 客观 分 
析 法 与 非 对 称 高 斯 函数 拟 合 ” ,将 8 d LST 数据 进 
行 重 构 , 为 了 与 NDVI 数据 时 间 分 辨 紊 统一 ,采用 平 
均值 合成 法 将 其 合成 为 16 d 数据 。 
1.2 地表 湿润 状况 估算 

Sandholt 等 指出 ,地 表 温 度 和 植被 覆盖 (NDVI) 
的 梯形 空间 可 以 反映 土壤 表层 的 水 分 含量 状况 ,所 
以 ,本 文采 用 TVDI 来 表征 地 表 湿 润 状况 ,TVDI 
基于 植被 指数 /地 表 温 度 的 梯形 特征 空间 计算 而 来 ， 
原理 如 下 : 


Ts — Ts ain (1) 
TS 一 
IP , Ts 为 每 个 像 元 的 LST; 四 为 对 应 像 元 湿 边 的 
LST 最 小 值 ;7s,, 为 对 应 像 元 湿 边 的 LST 最 大 值 
(图 1)。 
干 边 和 湿 边 的 拟 合 是 计算 TVDI 的 关键 环节 。 
首先 ,将 NDVI 的 值 域 0 ~1 划分 为 100 个 节点 , 节 


TVDI = 


点 之 间 相 差 0.01; 然 后 ,以 每 个 节点 为 中 心 左右 分 
别 扩展 0. 005 个 单位 形成 一 个 区 间 , 并 在 区 间 内 查 
询 LST 的 最 大 值 和 最 小 值 ;最 后 ,分 别 以 节点 值 序 
列 为 自 变 量 ,LST 的 最 大 值 和 最 小 值 序列 为 因 变 量 
进行 一 元 线性 回归 拟 合 ,从 而 得 到 干 边 和 湿 边 ,再 根 
据 式 (7 ) 计 算出 每 个 像 元 的 TVDI。 本 文中 干 湿 边 
拟 合 均 采 用 Theil-Sen 趋势 法 ,以 此 剔除 数据 中 存在 
的 异常 值 ,使 得 拟 合 的 干 湿 边 方程 更 为 稳健 。TVDI 
值 域 在 0 ~ 1 之 间 ,TVDI 越 大 ,地 表 越 干燥 ,TVDI 越 
小 ,地 表 越 湿润 。 


NDVI 


图 1 NDVI 与 LST( 地 表 湿 度 ) 构 建 TVDI 示意 图 
Fig.1 The TVDI configured by NDVI and LST 


根据 上 述 原 理 利用 16 d LST 和 NDVI 数据 计算 
得 到 16 d TVDI 数据 。 由 于 地 表 湿 润 状况 年 内 最 大 
值 不 能 反映 出 该 年 份 地 表 湿 润 的 总 体 状况 , 故 采用 
平均 值 合成 法 合成 TVDI 的 年 际 数 据 。 值 得 注意 的 
是 ,TVDI 的 计算 相当 于 是 对 LST 进行 最 大 限度 标准 
化 ,标准 化 的 上 下 限 不 是 一 个 定 值 ,而 是 随 着 NDVI 
的 不 同 发 生变 化 ,所 以 ,TVDI 与 NDVI 并 非 同 源 数 
据 , 可 以 反映 出 地 表 温 润 状况 与 植被 覆盖 之 间 的 关 
系 。 
1.3 Theil-Sen 趋势 

Theil-Sen 趋势 既 不 受 异 常 值 影响 ,也 不 用 服从 
一 定 的 分 布 , 结 果 更 为 科学 、 可 信 '*  。 所 以 本 文选 
用 这 种 方法 进行 趋势 分 析 , 并 采用 非 参 数 检验 法 
Mann-Kendall (M-K) "对 其 趋势 进行 检验 。 计 算 
公式 如 下 : 


w 


slope = Median ( 一 ) pya 


(2) 


S= > 5 sign (x; — x,) (3) 


i=lj=i+l 
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1, x,-x,>0 0.70 y=-0.004 5x+0.681 1 
P<0.05 
sign = (x, = x;) = 0, x,-x, =0 (4) 0.65 È 
-1, x -x; <0 = (ies SSS Sa Sea ie 
2 y=-0.006 3x+0.500 2 Z Ne 
人 B 0.55} P<0.001 AX i 
, S>0 E Ei 
VVar(9) 0.50 F —— TVDI —o— NDVI 
— TVDI slope NDVI slope 
Z = 0 9 S =0 (5) 0.45 1 1 1 1 1 1 1 1 1 1 1 L 1 
S41 2001 2003 2005 2007 2009 2011 2013 
= =0 年 份 
VVar(S) 
图 2 黄土 高 原 NDVI 与 TVDI 的 时 间 分 布 
P=2[1 (121)] (6) Fig.2 Temporal distribution of NDVI and TVDI 
E á in the Loess Plateau 
are 2 "7 ay (7) 覆盖 在 地 表 湿 润 状况 增加 的 影响 下 呈 增 加 趋势 。 年 
全 :一 一 t = y wy x 

. ra Bx TVDI 与 NDVI 时 间 序 列 的 Hurst 指数 分 别 为 


式 中 :slope 为 Theil-Sen 趋势 ;ww; 为 时 间 序 列 数 
据 ; 4,4; 为 时 间 序 列 对 应 的 年 份 ; 为 时 间 序列 的 时 
间 长 度 ;5 X M-K 相关 系数 ;Z 为 M-K 显著 性 统计 
值 ;P 为 显著 性 统计 值 对 应 的 概率 值 。 

本 文采 用 尸检 验 值 对 趋势 进行 检验 ,趋势 的 显 
著 性 被 划分 为 以 下 级 别 : 


R1 趋势 及 其 显著 性 划分 


Tab.1 Classification of the trend and significance 


趋势 P 检验 值 级 别 

slope >0 P<0.01 极 显著 增加 /改善 
0.01 <P <0. 05 显著 增加 /改善 

0.05 <P <0.1 弱 显 著 增 加 /改善 

P>0.1 不 显著 增加 /改善 

slope <0 P<0.01 极 显著 减 小 /退化 
0.01 <P<0.05 显著 减 小 /退化 

0.05<P<0.1 弱 显 著 减 小 /退化 

P>0.1 不 显著 减 小 /退化 

2 结果 与 分 析 


2.1 植被 覆盖 与 地 表 湿 润 状况 时 间 变 化 特征 
2001—2014 年 整个 黄土 高 原 NDVI 呈现 出 增加 
趋势 , 增 速 为 6.3% + (10a) “'(P <0.001) ,TVDI Œ 
现 出 下 降 趋 势 , 降 速 为 -4.5% + (10a) 
(P <0.05) ;年 平均 NDVI 均 小 于 0.6 ,年 平均 TVDI 
LIK F 0. 6 ,两 者 随时 间 变 化 趋 近 于 0. 6 ,但 未 曾 越 
过 这 个 值 (图 2), 且 两 者 在 时 间 变 化 上 呈 负 相关 , 相 
关系 数 为 -0.68(P<0.01) 。 表 明 , 黄 土 高 原 植被 


0.501 2 和 0.505 6 ,未 来 均 呈 微弱 的 持续 性 变化 , 即 
可 能 逐渐 接近 某 个 平衡 点 ,因此 ,0.6 可 能 是 NDVI 
和 TVDI 的 共同 平衡 点 ,两 者 越 接 近 这 个 值 ,黄土 高 
原 植被 覆盖 与 地 表 湿 润 状况 可 能 整体 进入 一 种 稳 
态 , 即 两 者 相互 依存 且 不 会 相互 耗 减 。 值 得 注意 的 
是 ,2012 年 以 后 ,植被 覆盖 远离 0. 6 这 个 平衡 点 , 植 
被 覆盖 出 现 减 小 趋势 ,从 而 导致 植被 保持 水 土 的 能 
力 减弱 ,地 表 湿 润 状况 可 能 会 出 现 减 小 态势 。 
2.2 植被 覆盖 与 地 表 湿 润 状况 空间 分 布 特征 

黄土 高 原 植被 覆盖 呈现 出 由 东南 向 西北 逐渐 减 
小 的 分 布 特征 ,高 值 区 集中 于 秦 巴 山地 区 和 燕山 一 
太行 山 山区 ,NDVI 均值 分 别 为 0.77 和 0. 73 ,该 区 
域 主要 分 布 着 落叶 与 常 绿 阔 叶 林 ,植被 茂密 ,其 次 是 
汾 渭 盆 地 和 农业 一 草原 过 渡 区 ,NDVI 分别 为 0.72 
和 0.56, 其 中 汾 渭 倪 地 属于 地 形 平坦 的 农业 区 ,分 
布 着 大 面积 农作物 ;农业 一 草原 过 渡 区 中 台 梁 山 以 
落叶 阔 叶 林 为 主 ,子午 岭 一 黄龙 山 属于 人 为 多 次 破 
坏 后 的 天 然 次 生 林 ,处 于 森林 草原 一 半 干 旱 草 原 过 
渡 地 带 , 六 盘山 系 典 型 草原 一 落叶 阔 叶 林 过 渡 区 ; 男 
外 ,江河 源 区 .河套 平原 .宁夏 平原 植被 覆盖 较 高 ,其 
中 江河 源 区 分 布 着 灌木 .草本 植被 .高 山 草 旬 ,河套 
平原 和 宁夏 平原 均 为 一 年 一 熟 粮食 作物 和 耐寒 经 济 
作物 。 低 值 区 集中 于 黄土 高 原 西北 部 的 非 农业 区 ， 
主要 为 草原 . 沙 地 和 荒漠 ,植被 稀 丽 ,NDVI 小 于 
0.45( 图 3 和 图 4) 。 

黄土 高 原 地 表 湿 润 状况 也 呈现 出 “东南 高 n 
北 低 ”的 分 布 特征 ,大 体 上 可 以 看 出 植被 覆盖 与 地 
表 湿 润 状况 具有 正 相 关 性 , 即 地 表 湿 润 状况 大 ,植被 
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(a) NDVI 


图 例 


EE < 0.3 
EE 0.3 ~ 0.4 
E 0.4~ 0.5 
æ 0.5 ~ 0.6 
0.6 ~ 0.7 
EE -0.7 


0 


100 200km 
— 


注 :生态 区 代码 I-10: 燕 山 - 太行山 山地 落叶 阔 叶 林 生态 区 ;I -11: 汾 渭 倪 地 农业 生态 区 ;1 -12: 黄 土 高 原 农 业 与 草原 生态 区 ;1-15: 秦 巴山 地 


落叶 与 常 绿 阔 叶 林 生态 区 ; 工 - 1: 内 蒙古 高 原 中 东部 典型 草原 生态 区 ; 工 -2 :内 蒙古 高 原 中 部 - 陇 中 荒漠 草原 生态 区 ; 工 -3 :内 蒙古 高 原 中 部 草 


原 化 荒漠 生态 区 ;IILE- 1: 祁 连 山 森 林 与 高 寒 草 原生 态 区 ;IIL-4: 江 河源 


-甘南 高 寒 草 甸 草原 生态 区 。 下 同 。 
图 3 ”黄土 高 原 NDVI 与 TVDI 的 空间 分 布 


Fig.3 Spatial distribution of NDVI and TVDI in the Loess Plateau 
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到 4 黄土 高 原 不 同 生 态 区 的 NDVI 与 TVDI 

Fig.4 The difference of NDVI and TVDI of different 


ecoregions in the Loess Plateau 


覆盖 高 ,反之 则 低 。 但 是 ,两 者 仍然 存在 不 一 致 的 地 
区 ,主要 为 汾 渭 盆 地 和 了 豫 西南 洛 河 两 侧 , 均 为 地 势 较 
为 平坦 的 农业 区 ,主要 分 布 着 一 年 两 熟 或 两 年 三 熟 
连作 粮食 作物 ,该 区 域 易 受 人 类 活动 影响 且 作物 常 
年 消耗 地 表 湿 润 状况 ,致使 地 表 湿 润 状况 较 低 。 
2.3 ”地 表 湿 润 状况 与 植被 覆盖 变化 趋势 

采用 Sen 趋势 计算 了 黄土 高 原 TVDI 和 NDVI 
的 变化 趋势 。 根 据 TVDI 变化 趋势 ,地 表 湿 润 状 况 
增加 和 减 小 趋势 的 面积 分 别 占 80.6% 和 19.4% ,地 
表 湿 润 状 况 整体 上 呈现 出 增加 趋势 。 地 表 湿 润 状况 
增加 区 中 ,研究 区 面积 的 51.72% 为 无 显著 增加 ,其 
次 是 显著 增加 \ 极 显著 增加 和 弱 显 车 增加 ,面积 分 别 
占 13.4% 8.8% Fil 6. 7% ;空间 上 ,地 表 湿润 状况 极 
著 、 显 车 与 弱 显 著 增 加 区 集中 分 布 于 山西 北部 的 
吕梁 山 恒山 五 台山 等 山地 、 陕 北 富 野 河 与 前 套 平 


和 


原 之 间 区 域 .陕西 段 黄 河 两 岸 、 陕 北 清润 河 与 延 河 之 
间 以 及 甘肃 段 渭河 南北 两 侧 , 且 呈现 出 以 极 显 著 区 
为 中 心 ,显著 性 逐渐 向 周围 减弱 (图 Sa) 。 地 表 湿 润 
状况 减 小 区 中 ,研究 区 面积 的 17.4% 无 显著 减 小 ， 
极 显 车 .显著 与 弱 显 著 减 小 区 面积 均 小 于 1.0% , 合 
计 2.0% ;空间 上 ,地 表 湿润 状况 减 小 区 主要 分 布 于 
研究 区 东南 部 陕西 渭河 一 洛 河 两 侧 河南 洛 河 两 侧 、 
山西 浊 漳 河 及 其 源 区 与 丹 河 两 侧 , 以 及 西北 部 鄂 尔 
多 斯 高 原 西部 一 黄河 以 东 的 荒漠 草原 区 .宁夏 平原 、 
昔 水 河 与 长 城 之 间 和 洗 水 谷地 (图 Sa) ,反映 出 地 表 
湿润 状况 减 小 区 主要 集中 于 河流 两 侧 且 毗连 城镇 建 
设 用 地 ,地表 湿 润 状况 的 减 小 与 河流 水 量 的 减少 , 导 
致 土壤 水 分 补给 不 足 、 人 类 活动 (人 口 增长 .土地 利 
用 变化 等 ) 密 切 相关 。 

根据 NDVI 变 化 趋势 ,黄土 高 原 植被 覆盖 改善 
和 退化 趋势 的 面积 分 别 占 85.9% 和 14. 1% ,植被 覆 
盖 整 体 上 呈 改 善 趋势 。 植 被 覆盖 增加 区 中 , 极 显 著 
改善 和 无 显著 改善 面积 相当 , 分 别 占 总 面积 的 
31.7% 和 32.2% ,其 次 是 显著 改善 和 弱 显著 改善 ， 
面积 分 别 占 16.3% AS. 6% ;空间 上 ,植被 覆盖 显著 
改善 区 集中 于 吕梁 山 以 西 一 子午 岭 和 黄龙 山 以 北 一 
白 于 山 以 东 一 毛乌素 沙 地 及 其 以 南 的 区 域 ,主要 为 
退耕 还 林 草 重点 建设 区 ,其 次 分 布 于 陇 中 高 原 及 
其 以 东 地 区 (图 5b)。 植 被 覆盖 退化 区 中 ,无 显著 退 
化 面积 占 11.2% ,显著 退化 区 面积 占 3. 0% ;空间 
上 ,植被 覆盖 退化 区 主要 分 布 于 农业 区 。 表明 植被 
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(a) TVDI 变 化 趋势 (b) NDVI 变 化 趋势 
图 例 图 例 

ME 极 显著 减 小 人 5 3 m 极 显著 退化 

E 显著 减 小 ME 显著 退化 

号 弱 显 著 减 小 弱 显 著 退 化 

也 无 显著 减 小 国 无 显著 退化 

ME 极 显著 增 大 国王 极 显著 改善 

ME GSK i 

了 弱 显 著 增 大 


上 无 显著 增 大 » 


Fe 


图 5 黄土 高 原 TVDI 与 NDVI 变化 趋势 空间 分 布 
Fig.5 The spatial distribution of TVDI and NDVI trends in the Loess Plateau 


(a) 相关 分 析 (b) 对 应 关系 
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图 6 黄土 高 原 地 表 湿润 状况 与 植被 覆盖 相关 分 析 


Fig.6 The correlation between surface humidity and vegetation cover in the Loess Plateau 


盖 变 化 基本 上 由 地 表 湿润 状况 变化 驱动 ,但 是 仍 。” 分 布 于 太行 山 一 吕梁 山 一 子午 岭 一 黄龙 山 等 山区 ， 
然 存 在 不 一 致 地 区 (图 5b)。 
2.4 植被 覆盖 与 地 表 湿 润 状况 的 时 空 关系 表 2 植被 覆盖 与 地 表 湿 润 状况 变化 趋势 类 型 划分 
2001 a, 年 黄土 高 原 TVDI 与 a 平均 相 Tab.2 Classification of the change trend of vegetation 
a =g IBJ JAN amn 
KABA - 0. 19, 且 负 相关 面积 占 研 究 区 的 Pre er See 
s on y pe nye IATA: 是 否 有 和 类 型 H 
73.9% WANA UME AES WM FAA REA VST gpi A EC 
现 出 正 相 关 ,通过 两 个 要 素 变化 趋势 的 空间 赦 加 ,得 ae i 改善 。 (图 6b) /% 
选 出 地 表 湿润 状况 变化 ST 与 植被 覆盖 变化 VT) RN eae 2 
、 k VT>ST RAAF 2 42.1 
对 应 关系 ( 表 2, 图 6b) ,从 而 反映 出 退耕 还 林 工 程 gx Ew 
以 来 两 者 的 变化 是 否 有 利于 生态 环境 改善 。 退化 增 思 ”一 - 较 有 利于 4 ”9.6 
5 
6 


cover and surface humidity 


(1) 植被 覆盖 改善 一 地 表 湿 润 状况 增加 且 其 趋 。 。” 退化 mN IVTI<ISTI BRAF 2.4 
SL VT < ST 的 区 域 (代码 1) 面积 占 28. 9% ,主要 ABA AE = 
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反映 出 14 a 间 植被 覆盖 与 地 表 湿 润 状况 具有 相互 
促进 作用 ,植被 覆盖 的 增加 提高 了 植物 水 源 涵养 
能 力 ,有 利于 生态 环境 改善 ; 而 植被 覆盖 改善 一 地 
表 湿 润 状况 增加 且 其 趋势 呈 VT > ST 的 区 域 (代码 
2) 达 42.1% ,分 布 于 退耕 还 林 重 点 区 5 ,反映 出 
该 区 域 植 物 生 长 速率 大 于 水 源 涵养 速率 ,短期 内 
生态 环境 处 于 改善 状态 ,长 此 以 往 可 能 会 导致 地 
表 水 分 耗 减 至 植物 生理 需 水 水 平 以 下 ,引起 生态 
再 次 退化 。 

(2) 植被 覆盖 改善 一 地 表 湿 润 状况 减 小 的 区 域 
(代码 3) 面 积 占 14.9% ,分 布 于 西北 部 鄂尔多斯 高 
原 西部 一 黄河 以 东 的 荒漠 草原 区 ,该 区 域 地 表 湿 润 
状况 基数 小 , 莹 发 量 大 ,加 之 植被 趋 于 改善 ,加 重地 
表 水 分 耗 减 ,生态 环境 不 易 在 人 类 活动 影响 下 改善 ， 
建议 采用 适当 灌溉 和 栽植 抗旱 植被 相 结合 的 方式 提 
高 生态 环境 质量 。 

(3) 植被 覆盖 退化 区 (代码 4,5,6) 面积 
14.1% ,主要 分 布 于 河套 平原 、 关 中 平原 、 汾 河谷 地 、 
涅 水 谷地 等 农耕 区 ,农耕 区 在 人 类 活动 影响 下 ,土地 
性 质 易 发 生变 化 ,引起 植被 覆盖 降低 ,如 粮食 作物 转 
换 为 经 济 作 物 (果园 ) ,植被 覆盖 的 退化 也 会 增加 土 
壤 侵 刨 ;而 地 表 湿 润 状况 却 在 不 合理 灌溉 措施 下 呈 
现 出 不 稳定 变化 趋势 。 该 区 域 属于 重要 的 粮食 生产 
基地 ,建议 恢复 原 有 作物 类 型 或 约束 作物 性 质 发 生 
变化 ,并 定量 核算 作物 灌溉 需 水 量 ,以 维持 作物 水 分 
供需 平衡 。 


3 讨论 


本 文采 用 TVDI 来 表示 地 表 湿 润 状况 大 小 , 通 
过 计算 地 表 湿 润 状况 和 植被 覆盖 年 内 变化 来 分 析 
TVDI, 反 映 地 表 湿 润 状况 的 可 靠 性 。 将 黄土 高 原 植 
被 覆盖 和 地 表 湿润 状况 年 内 数据 进行 Whittaker 平 
滑 和 传 里 叶 插 值 增 采 样 (356 d) ,一 定 程 度 上 剔除 了 
数据 噪声 ,提高 了 信 噪 比 (图 7) 。 

地 表 湿 润 状况 年 内 季节 变化 特征 与 张蕾 等 的 西 
北 地 区 土壤 水 分 变化 特征 ”高度 吻合 ,差异 性 主要 
表现 在 变化 特征 点 所 对 应 的 时 间 点 不 同 。 本 文 地 表 
湿润 状况 最 高 值 出 现在 一 年 内 的 1 月 上 旬 ,4 月 下 
旬 地 表 湿 润 状况 降 至 最 低 值 ,而 张 蓄 等 分 析 土 壤 水 
分 12 月 中 旬 至 1 月 较 高 ,5 月 下 旬 至 6 月 上 旬 降 至 
最 小 ,本 文 地 表 湿 润 状况 最 低 值 对 应 时 间 较 张 一 
提前 了 1 个 月 以 上 。 原 因 主 要 是 数据 的 时 间 分 辨 率 


0.6, 一 一 植被 覆盖 ”一 一 土壤 湿度 -0.5 
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图 7 黄土 高 原 植 被 覆盖 与 地 表 湿 润 状况 
的 年 内 季节 变化 特征 


Fig.7 The seasonal variation of vegetation cover and 


surface humidity in the Loess Plateau 


不 同 引 起 的 ,本 文采 用 16 d 数据 进行 分 析 ,而 张 董 
采用 旬 值 数据 ,分 析 中 体现 了 植物 生长 耗 水 引起 地 
表 湿 润 状况 变化 ,虽然 其 时 间 分 辩 率 较 高 ,但 缺乏 植 
被 覆盖 数据 作为 支撑 ,特征 时 间 节 点 与 实际 植物 生 
长 情况 不 符 。 从 植被 覆盖 年 内 季节 变化 特征 来 看 ， 
植物 生长 季 始 期 为 4 月 下 旬 ,与 地 表 湿 润 状况 最 低 
值 时 间 点 契合 度 很 好 (图 7) 。 韦 振 锋 等 分 析 得 出 陕 
甘 宁 黄土 高 原 大 多 数 植物 生长 季 始 期 为 3 一 4 
HES ,李强 等 认为 生长 季 始 期 为 3 月 下 旬 至 5 月 下 
°° ; 谢 宝 妮 等 认为 生长 季 始 期 为 4 月 上 旬 至 5 月 
FA ,证 明 TVDI 用 来 表示 地 表 湿 润 状况 变化 特 
征 具 有 很 高 的 可 信 性 。 

虽然 TVDI 具有 一 定 的 可 信 性 ,但 是 本 文 仍然 
存在 不 足 之 处 :利用 TVDI 反映 地 表 湿 润 状况 大 小 
未 经 实测 数据 的 时 空 验证 ,也 未 将 其 转换 为 土壤 水 
分 体积 分 数 , 后 期 数据 应 用 较为 困难 。 而 土壤 水 分 
时 空 变化 机 理 复 杂 , 以 及 传统 反 演 模型 具有 时 空 限 
制 .不 稳定 等 缺陷 。 在 此 基础 上 查阅 大 量 文献 ,发 现 
高 光谱 遥感 数据 为 土壤 水 分 反 演 提供 了 大 量 的 样本 
数据 , 且 反 射 波谱 曲线 的 细微 差别 可 以 反映 出 土壤 
水 分 的 微小 变化 ,使 得 单一 或 多 波段 遥感 数据 探测 
不 到 的 时 空 分 布 规律 ,在 高 光谱 遥感 中 得 以 实 
现 ”) 。 另 外 ,采用 机 器 学 习 方法 可 构建 任意 的 非 线 
性 关系 ,尤其 是 引起 广泛 关注 的 次 度 学 习 方 法 , EE 
感 影像 分 类 ,预测 与 反 演 等 方面 表现 出 强大 的 应 用 
潜力 。 因 此 ,今后 的 研究 将 侧重 于 土壤 水 分 的 
高 兴 谱 遥感 反 演 ,基于 高 光谱 遥感 数据 和 实测 土壤 
水 分 数据 ,利用 深度 学 习 方法 建立 遥感 数据 与 实测 
土壤 水 分 之 间 的 最 优 模型 ,从 而 提高 区 域 土壤 水 分 
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反 演 精度 。 
4 结论 


(1) 时 间 上 ,2001 一 2014 年 黄土 高 原 NDVI 增 
速 为 6. 3% + (10a)-!, TVDI 降 速 为 -4.5% ， 
(10a) ,表明 黄土 高 原 植被 覆盖 与 地 表 湿 润 状况 
均 呈 增加 趋势 。 空 间 上 ,植被 覆盖 与 地 表 湿 润 状况 
均 呈 “东南 高 西北 低 ” 的 分 布 特征 ,大 体 上 可 以 看 
出 两 者 具有 正 的 空间 相关 性 。 

(2) 从 TVDI 和 NDVI 的 变化 趋势 来 看 ,黄土 
高 原 植被 覆盖 与 地 表 湿 润 状况 整体 上 均 呈 现 出 增 
加 趋势 ;空间 分 布 上 反映 出 ,植被 覆盖 变化 基本 上 
致 ,但 仍然 存在 差异 

性 区 域 。 

(3) 从 相关 系数 来 看 ,黄土 高 原 植被 覆盖 与 地 
表 湿润 状况 呈现 出 正 相 关 ; 从 变化 趋势 对 应 关系 来 
看 ,14 a 间 黄土 高 原 东北 一 西南 一 线 山地 植被 覆盖 
与 地 表 湿 润 状况 变化 有 利于 生态 环境 改善 ;人 类 活 
动 强烈 的 平原 地 带 极 不 利于 生态 环境 改善 ;值得 注 
意 的 是 ,退耕 还 林 重 点 区 是 未 来 生态 环境 变化 监测 
和 管理 的 重点 区 域 ,以免 引 起 再 次 退化 。 
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Spatiotemporal Evolution of Vegetation Cover and Surface 
Humidity Since Implementing the Grain for Green Project 
in the Loess Plateau 


ZHANG Chong, WANG Jing, LEI Tian-wang, MA Ling, SONG Dian-xing 
(Shaanxi Key Laboratory of Disaster Monitoring and Mechanism Modeling , Baoji University of Arts and Sciences , 
Baoji 721013, Shaanxi, China) 


Abstract: Based on the MODIS data and Landsat data during the period from 2000 to 2014, the temperature veg- 
etation drought index was used to reflect the surface humidity, and the temporal and spatial changes of surface hu- 
midity and vegetation cover on the Loess Plateau since implementing the Grain for Green Project were analyzed. 
The findings suggested as follows: T) The growth rate of normalized vegetation index in the Loess Plateau during the 
period from 2001 to 2014 was 6.3% + (10a) ~', the temperature vegetation drought index was decreased by 
-4.5% + (10a) ~', the vegetation cover and surface humidity were increased, the vegetation coverage was high in 
the southeast but low in the northwest, and it could be generally seen that there was a positive correlation between 
the surface humidity and the vegetation coverage; (2) According to the change trend of TVDI and NDVI, the vegeta- 
tion cover and surface humidity were in an increase trend in the Loess Plateau. The spatial change of vegetation 
cover was in accordance with the variation of surface humidity perfectly, but there was a regional difference; (3) The 
correlative coefficient indicated that there was a positive correlation between the vegetation cover and the surface hu- 
midity in the Loess Plateau. From the perspective of the change trend, the changes of vegetation cover and surface 
humidity in the northeast-southwest zone of the Loess Plateau in 14 years were beneficial to the improvement of eco- 
logical environment. Intensive human activities in the plain are very unfavorable to the improvement of ecological 
environment. It is worth to pay attention to the key areas where the Grain for Green Project is implemented because 
they are critical for monitoring and managing the ecological environment change in the future so as to avoid degener- 
ating again. 

Key words: vegetation cover; surface humidity; spatiotemporal evolution; Grain for Green Project; Loess Plat- 
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